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We determine the complete NLL running of the static potential associated with the locally 

1/2 BPS Wilson loop in M = 4 supersymmetric Yang-Mills theory. We present results for the 

SU(N C ) singlet as well as for the adjoint configuration and arbitrary N c at weak coupling. In order 

to derive the respective anomalous dimensions we perform a two-loop calculation in the Jv = 4 

supersymmetric version of the effective field theory pNRQCD. In addition we confirm the recently 

ON obtained fixed-order result for the singlet static potential generated exclusively by ladder diagrams 

O 

£N) to the third order in the t'Hooft coupling. We also give an explicit expression for the logarithmic 
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contribution of all non-ladder diagrams at this order. 
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I. INTRODUCTION 



M = 4 supersymmetric Yang-Mills (SYM) theory in four spacetime dimensions is com- 
monly believed to be integrable in the large- N c (planar) limit. The reason is connected to 
the conjectured duality between certain pairs of conformal quantum field and string theory 
models, the so called AdS/CFT correspondence p~H3]. Notably it maps observables in planar 
M = 4 SYM theory to dual observables in type IIB superstring theory on an AdS^ x S 5 
background with vanishing string coupling. This string theory in turn belongs to a class of 
two dimensional (worldsheet) models, which appear to be integrable [I]. 

The key feature of Af = 4 SYM theory is its large symmetry content and in particular its 
conformal symmetry, which is necessary for the AdS / CFT duality to apply. As a consequence 
of the conformal invariance the beta function of the M = 4 SYM gauge coupling g vanishes 
to all orders. Correspondingly also the t'Hooft coupling A = N c g 2 , which is kept finite in 
the large- N c model, is scale independent. 

Despite the assumed integrability finding exact solutions in terms of algebraic functions 
of the coupling constant for observables like scattering amplitudes or expectation values 
of gauge invariant Wilson loops in planar = 4 SYM theory is far from trivial. A strong 
motivation to continue the effort in this direction is the hope that some (qualitative) features 
of such solutions are universal for a wider class of four- dimensional gauge theories including 
more realistic ones like QCD. Evidence for some universal behaviour has e.g. been found 
for gluon scattering amplitudes. Likewise in some instances M = 4 SYM results might 
possibly be regarded as the first approximation of the respective results in less symmetric 
Yang-Mills theories. On the other hand it is also conceivable to get new insigths into string 
theory/quantum gravity by studying exact results derived by means of field theory methods 
in M = 4 SYM theory due to the AdS / CFT correspondence. In order to formulate and test 
an all-orders ansatz for an observable in (planar) M = 4 SYM theory it is crucial to know 
the perturbative results at weak and strong coupling. The first can be computed within 
perturbative quantum field theory, the second within perturbative string theory exploiting 
the AdS/CFT duality. 

In this work we consider the energy of two static sources, which are in the fundamental 
(quark) and antifundamental (antiquark) representation of SU(N C ), respectively, and are 
separated by the spatial distance r. The system can either be in a SU(N C ) singlet or an 
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SU(N C ) adjoint multiplet configuration, the latter of which interacting with the Af = 4 SYM 
fields to form a physical quantity. In particular we focus on the singlet static energy given 
by 

E s (r)=\\m lln(Wb>, (1) 

1 — >oo ± 

where Wu is a rectangular representative of the class of 1/2 BPS Wilson loops 

W c = ^ Tr V exp ( - ig£ dr{A^ + $ n |£|)) (2) 

with the path C = □ spanned by the spacetime points x\ = (T/2,r/2), xi = (T/2, — r/2), 
yi = (-T/2,r/2) and y 2 = (-T/2, -r/2). The definition of the Wilson loop, Eq. @, is 
motivated in Ref. [5] with the static quarks modeled by W-bosons. We will often refer to 
the corresponding energy in Eq. ([T]) as the static potential. It is among the most prominent 
observables studied in the context of gauge/string duality. The main difference to the QCD 
analog is that the sources not only couple to gluons, but also to the six scalar fields $j of 
four-dimensional J\f = 4 SYM theory through the $ n = <fr • n term in Eq. ([2]), h being a 
six-dimensional unit vector (n 2 = 1). 

The perturbative calculation of the M = 4 SYM static potential in terms of quark- 
antiquark interaction diagrams relies on the (full theory) Lagrangian 1 

^ = C s tat + £-N=A ■ (3) 

Eq. ^ fixes the interaction terms of the static (quark) sources ip and x i n the fundamental 
and antifundamental representation, respectively, 

C s t a t = ^(ido - gA - g$ n )ip + x f (id + gA T Q - g$ T n )x , (4) 

and £v=4 denotes the usual Lagrangian of M = 4 SYM theory without external sources in 
four-dimensional Minkowski space [8] , 

C N=A = - \F; u F^ a + l - ^(D^TiD^y + h^D^.r + ■■■ . (5) 

1=1 

1 The analogy of the static potential to a quark-antiquark scattering process on the diagram level in mo- 
mentum space is actually only evident for a number of gauges including the generalized Lorentz gauge, 
cf. Ref. [S] for the QCD case. We will therefore restrict ourselves to the latter in the following discussion. 
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The tyj (j = 1,..,4) are four Majorana four-spinors living in the adjoint representation 
of SU(N C ) just like the gluon field N 1 and the scalar field $j (i = 1,..,6). For the 
gauge covariant derivative acting on these fields we adopt the following (sign) convention: 
{D^X) a = d^X a - gf abc A b ^X c , where X = A", 

A number of perturbative calculations of the static potential in (planar) M = 4 SYM 
theory in the weak as well as in the strong coupling regime can be found in the literature. 
For A ^> 1 present string theory calculations |9j [10] provide the first perturbative quantum 
correction of relative C(A -1 / 2 ) to the classical string result [H] See also Ref. [T3] for 
smooth generalizations of the rectangular Wilson loop to other shapes like the circular Wilson 
loop. A recent calculation of the M = 4 SYM static potential at weak coupling (A <C 1) 
considers diagrams of ladder type only [H]. Besides the explicit result up to (absolute) 
0(X 3 ) Ref. [13] also shows how a summation of all ladder diagrams can be achieved by 
solving a simple Schrodinger problem. While ladder diagrams are all there is at the tree, 
0(A), and one-loop level, C(A 2 ) [151 HE], diagrams with other topologies contribute at two 
loops, C(A 3 ), and beyond 2 . In Ref. [T7] a systematic generalization of the results of Ref. [H] 
in terms of an expansion around the ladder result for the static potential is discussed. As a 
first step the authors computed the leading logarithmic correction to the ladder limit. 

In the present work we are concerned with the weak coupling limit of the M = 4 SYM 
static potential. The strictly perturbative approach in terms of on-shell amplitudes derived 
from the Lagrangian in Eq. ^ is spoiled by infrared (IR) singularities. They occur first at 
one loop and do not cancel among each other in the individual terms of the loop expansion. 
A similar problem in the perturbation series of the QCD static potential has been known 
for a long time. In four spacetime dimensions the perturbative QCD corrections start to 
exhibit IR divergences at the three- loop [18] and in three dimensions at the two-loop level [6] . 
The reason for the divergent behaviour in the IR is that the strictly perturbative approach 
fails to allow for effects at the so-called "ultrasoft" scale ~ A/r. This scale parametrizes 
the difference between the binding energy of (intermediate) adjoint and singlet states of 
the quark-antiquark system. The ultrasoft scale does not appear in single quark- ant iquark 
interaction diagrams. It is rather generated by a resummation (exponentiation) of certain 
contributions of loop diagrams to all orders and acts as a physical IR regulator. 

2 Note that in contrast to Ref. [14] we stick to the standard counting of loops in QCD (scattering) diagrams. 
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In Ref. [H] a corresponding resummation is explicitly performed for ladder diagrams and 
finite results through next-to-next-to-leading order (NNLO), i.e. 0(X 3 ), are obtained. These 
results are now non-analytic in the coupling constant A since they contain logarithms of the 
ratio between the soft scale ~ 1/r and the ultrasoft scale ~ A/r. 3 They lead to the typical 
In A terms, which are strongly connected to the logarithmic IR singularities of the strictly 
perturbative calculation. 

A more systematic solution to the problem of IR singularities in the QCD static potential 
has been proposed in Ref. [TH] and has been successfully applied not only in four, but 
also in three spacetime dimensions [20] • It is based on the effective field theory (EFT) 
potential nonrelativistic QCD (pNRQCD) [21] (for a review see Ref. [22] ). This framework 
explicitly contains degrees of freedom living at the ultrasoft scale, most prominently ultrasoft 
gluons. It allows to systematically calculate the ultrasoft contributions missing in the strictly 
perturbative calculation from a finite number of EFT diagrams at every order in the coupling 
constant. Most notably the ultraviolet divergences of the pNRQCD results exactly cancel 
the IR divergences from the respective full theory (QCD) diagrams, which (in the static 
limit) only involve the soft scale 1/r. 

The analogous EFT for M = 4 SYM theory was introduced by Pineda in Ref. [23]. Also 
the complete next-to-leading order (NLO), i.e. C(A 2 ), result of the singlet static energy E s {r) 
at weak coupling was derived there for the first time. The result also holds for arbitrary iV c . 
In addition, Ref. [23] makes use of the renormalization group (RG) formalism to determine 
the leading logarithmic (LL) running of the static potential from the leading ultrasoft UV 
divergence in the planar limit. To some extent this already represents an all-orders result 
and demonstrates the power of the EFT approach. 

It is the aim of the present work to compute the next-to-leading logarithmic (NLL) RG 
evolution of the static potential in this EFT framework for arbitrary N c . With the result we 
will be able to deduce also the full NNLO expression for the singlet static potential in the 
planar limit, including the contributions from non-ladder diagrams, up to a (soft) constant 
and thus generalize the results of Refs. [HI ITT] . 

The outline of this paper is as follows. In Sec. [TTJwe briefly review the theoretical frame- 

3 The terminology for the scales (soft, ultrasoft) actually comes from the more physical situation of dynami- 
cal quarks with a mass, which is large but finite and defines a third so called "hard" scale m ^> 1/r ^> A/r. 
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work we are working in and explain the ultrasoft two-loop calculation we have performed. 
Explicit results for individual diagrams are given in App. [B] In Sec. Ill we solve the renor- 



malization group equations (RGE's) derived in Sec. |TT] allowing for the NLO soft matching 
conditions for the singlet and adjoint static potentials. We then give the complete NLL re- 



sults for the respective static energies. Section IV is concerned with the large- N c limit of our 



results. After computing also the soft two-loop matching condition from ladder graphs we 
finally present the full large- N c NNLO result for the singlet static energy up to an unknown 
constant coming from soft two-loop non-ladder diagrams. We conclude in Sec. |Vj 

II. THE EFT CALCULATION 

The EFT we are employing in this work was introduced in Ref. [22] and is based on 
pNRQCD. It exploits the hierarchy between the soft and the ultrasoft scale in the weak 
coupling limit, £' s ~A/r^l/r and makes the factorization between these scales manifest. 
This is achieved by the "multipole" expansion of the effective Lagrangian in powers of 
the distance r. Fluctuations with energies and/or momenta ~ 1/r are integrated out. The 
remaining degrees of freedom are massless gluons, scalars and fermions with ultrasoft energies 
and momenta and the quark-antiquark pair with ultrasoft energy and zero momentum. The 
latter is conveniently projected on the physical singlet (S) and adjoint (O) field. In analogy 
to QCD we will often refer to the SU(N C ) adjoint as the "octet" in the following. 

At leading order in the multipole expansion the singlet-octet sector of the EFT La- 
grangian [23] reads 

£ us = Trjst (id - V s (r)) S + 0* (iD - V (r)) O j 

- 2gV*(r)Tr {S^O + S^O*} - gV* (r)Tt {O* O}} + 0{r) . (6) 

S and O are operators in color space and related to the normalized singlet and octet fields 
by 

1 T a 

' S, = ^O a , (7) 



where the T a are the generators of the fundamental representation of SU(N C ) (Tp = 1/2). 
The EFT Lagrangian of course respects all symmetries of the full theory Lagrangian, Eq. pi). 
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The terms in Eq. (|5]) retain their form for the ultrasoft gluons, scalars and fermions. Accord- 
ingly the Feynman rules for the interaction of these fields among each other are the same 
as in full M = 4 SYM theory. The EFT Feynman rules from the Lagrangian in Eq. ^ are 
listed in Appendix |A} 

The Wilson coefficients Vi depend on the distance r and on the renormalization scale v. 
They are typically fixed at the soft scale [v ~ 1/r) by matching to the full theory. We will 
refer to the coefficients 

V s {r) = -2C F — — , V {r) =21 — Cp J — — , (8) 

as the (soft) singlet and octet static potentials, respectively {Ca = N c , Cf = (N^ — l) / (2N C )). 

2 

Tree level matching gives oty s = «y o = a = |-, V$ = V$ = 1. The potential V s / Q is to be 
distinguished from the observable static energy 4 

E s/0 (r) = V s/0 + 5E:j . (9) 

The Wilson loop definition of the static octet energy differs from the singlet energy, Eq. (JT|), 
by an insertion of the SU(N C ) generator T a on each of the spatial lines of the rectangular 
loop, see e.g. Ref. [22] for the QCD analog. 5E US is the ultrasoft contribution to the static 
energy. We will compute it below for the singlet as well as for the octet at the two-loop 
level. Eq. ^ holds for renormalized as well as for bare quantities on the right hand side. 

In calculations of SE US the ultrasoft scale is represented by AV = V Q — V s . In fact the 
potential coefficients V s and V Q only appear in this combination. At NLO in the multipole 
expansion the operators of £us are proportional to the the distance vector r and come with 
at least one power of the coupling constant g, see Ref. [23] • Since rotational invariance 
requires an even number of NLO operators their contribution to the static energy scales like 
aAV 3 r 2 ~ a 4 /r. The NNLO operators (oc ar 2 ) start to contribute at the same parametric 
order. This is beyond the precision we are aiming for in the present work. We can therefore 
safely neglect operators of dimension greater than four and work only with the LO operators 
in Eq. @. 

4 As mentioned above we will sometimes also use the term "potential" for the static energy. We will however 
be more precise, whenever the distinction between the Wilson coefficient V s / a and the observable energy 
E s / is important. 
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Let —iYj s / (E) denote the sum of all singlet /octet 1PI self-energy diagrams in the EFT. 
Then the full propagator of the interacting singlet /octet field takes the form 



% E^E s/o (r)^ iZ s / 

E - V s y™ - Z 8/o (E) * E - E s/o 



> i^hzs- (10) 



The arrow indicates the physical on-shell limit and \JZ S / Q is the singlet/octet wave function 
renormalization constant. In the following quantities without the "bare" label are under- 
stood to be renormalized. From Eqs. (|9| and (10) one finds 



8E U S J = Z s/o (6V s/o + H s/o {E = V s/o )) + 0(a 3 ) , (11) 
Z s/o = (1 - T,' s/0 {E = V s , )Y l + 0{a 2 ) (12) 

to second order in the loop expansion. The term T,' S ^ (E = V S ) denotes the derivative of the 
self-energy £ s / with respect to the external energy E evaluated at E = V s / and SV s / a = 
Vgfo £ ~ Vs/o is the counterterm of the singlet/octet potential in a given renormalization 
scheme. 

Before we turn to the actual calculation of 8E^j a we first consider possible higher order 
corrections to the Wilson coefficients V$ and V$ Q as they can mix into the result for the static 
energy. The coupling constant g does not receive corrections due to conformal invariance 
inherited from the full theory. 5 It has been stated already in Ref. [23] that there are no 
corrections to and V^ Q at 0(a). Indeed we have found vanishing LL anomalous dimen- 
sions by direct computation in the EFT. Inspecting the soft one-loop diagrams relevant for 
the 0(a 2 ) matching of and V^ Q one finds moreover that they possibly contain IR linear 
but no logarithmic divergences. As a consequence also the NLL anomalous dimensions must 
vanish. Summarizing we have (n = 0, 1, 2, . . . ) 

V* = l + 0(a 2 (a\na) n ) , V+ = 1 + 0(a 2 (a\na) n ) . (13) 



Thus mixing corrections to Eq. (11) are higher order and we will consistently use = 
V$ = 1 in the following. 

Throughout this paper we will work in Feynman gauge using dimensional regulariza- 
tion (d = A — 2e) and the MS renormalization scheme, unless explicitly stated otherwise. 



Wc have checked this statement explicitly at one loop level in the EFT. 
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FIG. 1: One-loop self-energy diagrams for the singlet (a) and octet (b) field. The dotted line 
represents the scalar the double line the octet and the bold single line the singlet. The crosses 
symbolize the O(r ) singlet-octet-scalar interaction from Eq. [fij The corresponding Feynman rules 
are found in Appendix \A\ 

We parametrize our results for the ultrasoft contributions to the static energies and the 
anomalous dimensions of the potentials as follows 

2 

5E™ = -C F AV ( - + — + 0{o?) ) + 0{r 2 ) , (14) 
SET = - C F ) AV ( ^ 6 W + g 6 00 + 0(a ») ) + 0(r 2 ) , (15) 

vj- V s = -C F AV ( - c?) + ^ c?) + G(a 3 ) ) + 0(r 2 ) , (16) 
^ K = - Cfe) AV ( £ c« + £ <f + 0(a 3 ) ) + 0(r 2 ) . (17) 



The latter are determined from the potential counterterms in Eq. (11) in the standard way. 



using vdajdv = —lea. Note that for large N c the octet expressions in Eqs. d8J), (15) 



and (17) are suppressed by two powers in the 1/N C expansion as compared to the singlet. 

At one-loop level the relevant singlet and octet self-energy diagrams are depicted in 
Fig. [TJ There are also two one-loop octet self-energy diagrams without an intermediate 
singlet propagator, which vanish however for E = V Q and are not shown. 6 The singlet 
diagram in Fig. [T^, has already been computed in Ref . [23] . Together with the octet diagram 
in Fig. [T]d we find 

6« = 41n( 2 ^)-4, 6 W = 4 ]b(=^)-4, (19) 
c?) = c« = 4 . (20) 

6 It is also possible to perform the calculation with small, but hnite E — V a as a regulator for intermediate 
IR divergences, cf. Ref. [21] ■ We have checked that the additional terms induced by this regulator cancel 
in our final results at one and two loops. 
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The relevant two-loop diagrams and their individual results after integration in d dimen- 
sions are shown in Appendix [Bj We have solved the integrals using integration by parts and 
partial fraction techniques. Expanding the respective self-energy expressions in e = (4 — d)/2 
and absorbing the UV divergences into the potential counterterms SV S / Q yields 



b^=4C A \n 2 



/2AV 



\ v 



+ ,4 + £)c^: 



(2AV 



4vr 2 

(12-8C(3))CU-— C F 



(21) 



6f=4C A ln 2 (^) + ((l6 



in 



+ (8C(3) - 36 + 4tt 2 )Ca + (64 - ^- )C, 



., y C A -32C F )ln(— — ) - 

47T 5 



.(2) 



.(2) 



4tt 2 - 

4 + -)r, 



16 



in' 
~3~ 



C^4 — 32 C p . 



(22) 
(23) 
(24) 



In QCD, in contrast, nonzero ultrasoft contributions to the singlet and octet static ener- 
gies and potential anomalous dimensions first occur at NNLO in the multipole expansion, 
i.e. at 0(r 2 ). Explicit results of state-of-the-art pNRQCD two-loop calculations at this 
order can be found in Ref . [24J . They are potentially relevant for calculations of the M = 4 
SYM static energy at 0(a 5 /r) or the "ordinary" Wilson loop without couplings of the static 
quarks to the scalar field. 



III. RG EVOLUTION AND MATCHING 

In order to perform the NLL resummation for the static potentials we have to solve the 



RGE's in Eqs. (16) and (17). For finite N c , i.e. nonvanishing octet potential, the Wilson 
coefficients V s and V Q mix into each other. Their NLL RGE's are coupled due to the linear 
dependence on AV = V a — V s . Therefore we better rewrite them in matrix form: 



d 

l 

dv 



"Vs" 


= M 


'V s ' 


te := 


'-A 


A' 


-Vo- 




-Vo- 




.-B 


B. 



(25) 



with the coefficients 



A = -C F (4- + U3 + n 2 )C A - 2 ) 

c A 



B 



-Ci 



a 



7T 



(( 6 



12C F ) 



a 



(26) 
(27) 
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The solution to Eq. (25) is given by 



(28) 



w l r/ifi t w \ Vs[llr) ~ 

= exp \M mvru) 
V (u)\ - W (l/r). 

where V s / {l/r) = V s / (r; u= 1/r) is understood. Carrying out the matrix exponentiation in 



Eq. (28) we obtain 
[V s 

[Vo 



[V 



1 



B-Al 
1 



BV 8 (l/r) - AV (l/r) + AAV(l/r) (rv) B - A 
B V.(l/r) - A V a (l/r) + B AV(l/r) {rv) B - A 



(29) 
(30) 



B — A 

What remains to be computed for the complete NLL solutions are the NLO matching 
coefficients V S)0 (l/r). The perturbative TV = 4 SYM theory calculation of the tree-level and 
the one-loop (crossed) ladder diagrams yields the bare potentials 



V S , B 

v 0jB 



2C F g 2 r 2e 1 T[l/2 - e] 



47T r 7T 1 / 2 e 

(C F - C A /2) 



C A g 2 r 



r[l/2-2e]r 2 [-e] 
87r 2 - e r[-2e]r[l/2-e] 



+ 0(g 4 



V s ,b + 0(g b ) . 



(31) 
(32) 



To extract the MS renormalized expressions we replace g — > gv\ exp[e (7^ — ln(47r))/2] in 



Eqs. (31) and (32), subtract the respective potential one-loop counterterms and take the 



e — > limit. In this procedure it is important to keep the 0(e) term of the LO MS expression 
for AV, which is inserted in the counterterms 8V s / . Otherwise we would miss a finite piece 
that must be subtracted from the bare potentials. Finally we set v s — 1/r and find 

2C F a 



V s (l/r) 



(l + 2C AlE - + 0(a 2 

\ 7T 



V.(l/r)= (C ^ 2Cf) + 20^ + 0(a»: 
r \ 7r 



(33) 
(34) 



Adding the ultrasoft one-loop terms of Eqs. ( 14 ) and ( 15 ) to Eqs. ( 29 ) and ( 30 ) , respectively, 



we can now write the static singlet and octet energies with NLL resummation as 

, L 4CV a 



[V s \ 

[Vol 



v 



[AV] LL (ln(2[AV] L V- 1 )-l), (35) 
NLL + 2(C A -2C F )a [Ay] LL ^ ln ( _ 2 [Ay] u.^ _ ^ (36) 



These expressions have a residual v dependence, since the ultrasoft logarithmic terms at 
0(a 3 ) or higher, e.g. the ones in Eqs. (21) and (22), are not included. In order to render 



the higher order ultrasoft logarithms small we should set v ~ in Eqs. (35) and (36). 
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Upon expansion in a Eqs. (35) and (36) yield the corresponding NLO expressions, 



7T 



(37) 



[E (r)} NL ° = ^ A ~ 2CF)a ( 1 + ^ [ H -2C A a) + lE - 1] ) , (38) 



r \ 7r 

where the v dependence of the soft and ultrasoft contributions has canceled completely. 



Moreover we can extract all LL (oc (a\na) n ) and NLL (oc a(a\na) n ) terms from Eqs. (35) 
and (36) by choosing v ~ AV ~ C A a/r. In particular we can determine the NNLO 



corrections to the static energies up to a real constant c^ LO : 



5[E S ( 
S[E ( 



AC 2 A C F a 3 
7r 2 r 



7T 



ln 2 {2C A a) + 1 + 2j E + — ln(2C» + c 



2Ci(C A -2C F )a 3 
7r 2 r 



x 



x 



ln^-20)+ 4 + 2 7E - 



8C f 



3 / 



2tt s 



(39) 
(40) 



ln(-2C A a) + C 



We find the same results if we add also the two-loop ultrasoft contributions in Eqs. (14) 



and (15) to Eqs. (35) and (36), respectively, replace [AV] by 

[AU] NL0 



C A a ( 2C A a f ,\ 



(41) 



and expand in a. To fix the constants c^ LO only the presently unknown 0(a 3 ) terms in 



the (scheme dependent) soft matching conditions Eqs. (33) and Eqs. (34) are missing. The 



respective ultrasoft components are already included in Eqs. (14) and (15). 



Obviously the expressions for the octet static energy in Eqs. (36), (38) and (40) have a 



nonzero imaginary part. It is related to the decay of the octet state into the singlet and 
massless quanta in the adjoint representation (gluons, scalars, Majorana fermions), 



lm[E {r)\ = -r o /2 



(42) 



whereas the real part is the physical energy. Since for v = AV the potential coefficient 
V Q (and AV itself) is real, we can determine the total decay width T directly from the 



ultrasoft contribution in Eq. (15). Note that every AV in Eq. (15) is accompanied by a 



positive imaginary infinitesimal due to causality: AV — > AV + ie, cf. the propagator (A4). 
We thus obtain 

4a 2 /. 2a 



( 1 + ^ [3C A \n(2aC A ) + (6 + 3 lE - vr 2 ) C A - 12C F ] + 0(a 2 



(43) 
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where we have used that Ca{Ca — 2CV) = 1. To put the energy and decay width of the octet 
quark- ant iquark system into a physical context one actually has to allow for the (ultrasoft) 
interaction with some gluon and/or scalar background needed to form an overall color singlet 
state, cf. Ref. [21] for the QCD case. This is however outside the scope of this paper. 



IV. THE LARGE-iV c LIMIT 

This section addresses the fixed order calculation of the static singlet energy in the large- 



N c limit. Eq. (14) already contains the ultrasoft terms through 0(X ). For the complete 



NNLO static energy we however also need the 0(X 3 ) term in the soft matching coefficient, 



Eq. (33). This requires a soft two- loop calculation. As a first step we have explicitly 
computed the two-loop ladder diagrams in perturbative M = 4 SYM theory. 7 We have 
evaluated the relevant momentum-space two-loop integrals using the techniques described 
in Ref. [6] and afterwards performed the Fourier transformation to (^-dimensional position- 



space. From the result we obtain the two- loop ladder contribution to be added to Eq. (31). 
In total we find 

yba re = _ A r* i r(l/2 - e) f £[1/2 - 2e]r 2 [-e] 

v s,lad An r nl/2 . t | * 87r 2 ^r[-2e]r[l/2-e] 

AV< 27rr( ~ 3e ~ 1)F3(1/2 - e) + r(1/2 - 3e)r3( ~ e) + 0(\*) \ (44) 
~ Xr 967r^(2e + l)r(-3e-l)r(l/2-e) + ° {X >J ' {U) 

The octet static potential vanishes in the planar limit so that AV = — V s . 

Since in our analysis we will treat ladder and non-ladder contributions separately, we 
also have to split the ultrasoft contributions and counterterms into two parts. This is easily 
possible, because full theory diagrams with a scalar or gluon (A ) propagator connecting 
two points on the same static quark line cancel in the sum. Accordingly only the ultrasoft 
diagrams with direct interaction among the scalars, gluons and/or Majorana fermions, i.e. 



7 In order to avoid pinch singularities we have only computed the maximally non-Abelian term oc C\Cf of 
the non-planar box (crossed ladder) diagrams. Although the total color factor of the non-planar diagrams 
vanishes in the large-iV c limit only the maximally non-Abelian term survives in the sum of all two-loop 
ladder diagrams after subtraction of the Coulomb pinches due to exponentiation. 
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diagrams (B7) and (B8), are of non-ladder type. Thus we have 

A 2 " 



$ E slad 



X T? us 

s,nonlad 



AV 



+ 0(r 2 



-AV 



A 

2t? 



X 2 



1 



7T 



3 V v 



+ 



7T 



57T 1 



to hr) +ln hr-)- 8 



+ 0(A 3 



+ 2C(3) 



+ 0(A 3 ) ^ + C>(r 



The associated counterterms read 

AV 

OV s j ad ■ 
^^s,nonlad ' 



4tt 2 
AV 



A 



+ 



A ", • o(x 3 



?7r 2 e 



+ 0(r 2 ) 



A 2 



+ C(A 3 ) +C(r 



(45) 
(46) 

(47) 
(48) 



To convert the bare expression Eq. (44) to the MS scheme we introduce the renormal- 



ization/factorization scale v by replacing A — > Xv 2e exp[e (7^ — ln(47r))] and subtract the 
counterterm 5V S i a d- Note that now subleading terms of AV^ not only in e, but also in A must 



be taken into account in Eq. (47). Before we can write down the NNLO result we therefore 



have to determine the NLO expression for V s = — AV, because the latter in turn is required 



in Eq. (47) as an input to reach NNLO precision. The outcome of this procedure is 

r2" 



V s i ad (v) 



A 



4-7T r 



1 



2^ 
2^ 



ln(e 7E rz/) + 



A 



■57T 



In [e tE rv] + \n(e lE ri/) 



w 
12 



+ 0(X J 



(49) 



Adding the ultrasoft ladder contribution of Eq. (45), where again the NLO expression for 



AV^ related to the first two terms of Eq. (49), must be inserted, we arrive at 

~ A 2 



A J A 

E St lad - - -. \ 1 + — — 2 

4:71 r Z7T Z 



+ C(A S 



In 



e 7£ A 
2tt 



1 



In 2 1 — — ) + In 



2tt 



2tt 



7T 

12 



(50) 



This expression exactly agrees with the ladder result for the static energy in Ref. 



As for the non-ladder part of the static energy only the analog to Eq. (49) is missing 



Its logarithmic term is however uniquely determined by requiring that the //-dependence 



exactly cancels the one of the ultrasoft non-ladder part in Eq. (46). This is equivalent to 



the postulation that the IR divergence in the soft non-ladder calculation must match the 



ultrasoft (1/e) UV divergence in Eq. (48). We can thus write 

■2 



V s nonlad\L / ) 



a r a 



47T r 87T 4 



7T 



ln(e 7B rz/) + c 



nonlad 



+ e>(A 2 



(51) 
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and together with Eq. (46) 

E s nonlad ~ 



A 



A 2 fvr 2 



An r 87r 4 



-ln(^ 
3 V 2tt 



+ c 



NNLO,soft 
nonlad 



7T 



- y + 2C(3) 



+ £>(A 3 



(52) 



We note that due to the additional factor 7r 2 /3 the single non-ladder logarithm at 0(X 3 ) 
is numerically dominant compared to the corresponding ladder logarithm. The sum of the 



ladder and non-ladder logarithmic terms oc A 3 (ln A)" in Eqs. (50) and (52) consistently agrees 



with the large- limit of Eq. (39 ). The real constant c^ z °'J oft can be obtained from the soft 
two-loop matching condition for the non-ladder potential V Synon i a d in the MS scheme, which 
is unknown at present. 

A corresponding soft two-loop calculation was performed in Ref. [14] . In contrast to our 
approach, however, the potential in Ref. [14J was computed on the S 3 sphere, where the 
static quarks are separated by an angle 5. The limit 5 — > then yields the static potential 
in flat space. In this limit 5 plays the role of the distance r and at the same time acts as IR 
cutoff. The (soft) results of Ref. [14] . 



V s 

' s,lad 



V s 

y s,nonlad 
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x , f 2s \ 
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57T 
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+ ln z 
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25 ^ 
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7T 
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+ C(A 3 



A A 2 Ttt 2 



1 3 4 



+ 0(A 3 



(53) 
(54) 



47r5 87r 4 

are therefore expressed in a (cutoff) scheme different from our MS scheme. Note that 



Eq. (53) contains a power- like IR divergence of the form l/<5 2 . It is canceled exactly by the 
corresponding ultrasoft expression in the 5-scheme [H] and does not have an equivalent in 



Eq. (49), because the latter has been derived using dimensional regularization. 



It is tempting to try to translate Eq. (54) to the MS scheme in order to fix the constant 
c mm°ad 0tt i OY alternatively translate Eq. (46) to the ^-scheme, by devising a dictionary for 
the logarithmic terms (in addition to the simple replacement rule 1/5 <H- 1/r) based on 
the knowledge of the respective ladder results in both schemes. We however refrain from 
doing so since the conversion of multi-loop results between different regularization schemes, 
in particular cutoff schemes and dimensional regularization, is a delicate issue and naive 
replacement rules can be misleading. Given the involved loop structure in our case a reliable 
translation prescription would require a dedicated study which goes beyond the scope of this 
paper. On the other hand a direct (MS) computation of c^^ oft from soft non-ladder two- 
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loop quark-antiquark interaction diagrams in flat space is extensive, but certainly feasible 
with present Feynman diagram technology. This work is left for the future. 
Summarizing the fixed order results in the large- N c limit we have 
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(55) 



Here we have parametrized the currently unknown terms by coefficients cf and explicitly 
expanded out the LL and NLL large- N c logarithms at 0(X A ) and 0(X 5 ) from Eq. (35) by 



choosing v ~ A/(47rr) for illustration. All the 7r 2 -enhanced NLL logarithms in Eq. (55) 
come from non-ladder diagrams, all other LL and NLL logarithms are generated by ladder 
diagrams. 



The one- loop logarithm in the first line of Eq. ( 55 ) was first determined by Erickson et al. 



in Ref. [15], the one-loop constant and all LL logarithms by Pineda in Ref. [23J. The 0(X 3 ) 



ladder result, Eq. (50), was first calculated by Correa et al. in Ref. |14] . Bykov and Zarembo 
computed in Ref. [17] the leading logarithmic correction to the ladder result, which agrees 



with our 0(X 3 ) non-ladder logarithm in Eq. (52). They however considered their result as 
incomplete and assumed that there are also other terms oc A 3 In A at higher orders in their 
expansion about the ladder result, which they had not computed. We have shown that this 



is not the case and Eq. (55) contains the complete LL and NLL, ladder as well as non-ladder, 



contributions. The remaining terms displayed in Eq. (55) are also new 



As already mentioned the NNLO constant, c^^°^ oft , can be computed in terms of soft, 
i.e. full theory, two-loop non-ladder diagrams. The NNLL coefficients cf NLL require the 
UV-divergences of the ultrasoft three- loop self-energy at 0(r°) and the ultrasoft one-loop 
self-energy at 0(r 2 ) [23], whereas for the determination of c n3lo the respective ultrasoft finite 
parts in addition to the soft three-loop calculation are needed. Finally for the complete 
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N 4 LO static energy four-loop soft and ultrasoft calculations have to be performed, where 
the ultrasoft four-loop O(r ) and two-loop 0(r 2 ) UV-divergences suffice to fix the N 3 LL 
coefficients cf 3LL . 



V. CONCLUSION 

The M = 4 SYM potential (energy) between two static sources in the fundamental 
representation of SU (N c ) is composed of soft and ultrasoft contributions. Using the effective 
theory proposed in Ref. [23] we have calculated the respective ultrasoft contributions for the 
singlet and adjoint color state at two loops. 

The main results of the paper are the NLL expressions for the M = 4 SYM singlet and 



adjoint static energy in Eqs. (35) and (36) as well as the large- N c fixed order expression for 
the singlet static energy in Eq. (55). For the renormalization scale choice v ~ C^a/V the 
former two expressions contain the complete resummation of logarithmic terms oc (alna) n 
and oc a(alna) n for arbitrary N c . The latter expression includes the contributions from 
all relevant diagrams to 0(X 3 ) except for an unknown constant. This constant can be 
determined from a purely perturbative (soft matching) computation of the two-loop non- 
ladder static quark- ant iquark interaction diagrams in flat space. We leave this calculation 
for the future. 
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Appendix A: EFT Feynman Rules 

This appendix lists the momentum space Feynman rules associated with the O(r ) EFT 
Lagrangian in Eq. (JH]). 

1. Singlet/octet vertices: 



■ o. 



(Al) 



6 1 C : -igV* d abc 



(A2) 



: gf 



abc 



(A3) 



2. Singlet/octet propagators: 

: if-v. + u (A4) 

= k"-K+ it (A5) 



Conventionally a left-right flow of the ultrasoft four-momentum fc M is understood. The 
dotted line represents the scalar the double line the octet, the bold single line the singlet 
and the spiral line the zero component of the gluon field (A ). The letters a, b, c denote color 
indices in the adjoint representation. The f abc and d abc are the totally antisymmetric and 
symmetric structure constants of SU(N C ), respectively. 
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Appendix B: EFT Diagrams 



In the following amplitudes we set = V$ = 1 and T F = 1/2. The expression for 



the scalar one-loop self-energy subgraph in the diagrams (B8) and (B16) can be found in 

Ref. ma. 



1. Two-loop singlet self-energy diagrams: 
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2. Two-loop octet self-energy diagrams: 



(Left-right mirror graphs are included in the expressions below.) 
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